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Abstract. We obtain a logarithmically sharp estimate for the space-analyticity radius of the solutions 
of the damped-driven 2D Navier— Stokes equations with periodic boundary conditions and relate this 
to the small scales in this system. This system is inspired by the Stommel— Charney barotropic ocean 
circulation model. 

Key words: Analyticity, Gevrey regularity, Navier— Stokes equations, dissipative length scales, 
Stommel-Charney model. 

AMS subject classification: 35B41, 35Q30, 37L30. 

1. Introduction 

It was shown in [15] (see also [3] , [12] ) that the solutions of the 2D Navier-Stokes equa- 
tions with periodic boundary conditions belong to the Gevrey class of analytic functions 
(if the forcing term does). Using the Gevrey regularity approach the following estimate 
for the spatial analyticity radius for the solutions that lie on the global attractor (or are 
near it) was obtained 

c\Q\ 1/2 

la ~ G 2 log G ' (L1) 

where G = ||/||l 2 I^I/ 1/2 i s the Grashof number and |0| = L 2 / r y is the area of the periodic 
domain ft = [0, L/7] x [0, L], 7 < I. 

Therefore, the Fourier coefficients Uk are exponentially small for \k\ ^> L/l a , and l a 
naturally forms a lower bound for the small dissipative length scale for the system (see, 
for instance, [EE]). 

There are other ways of estimating the dissipative small length scale for the Navier- 
Stokes system, for instance, in terms of the dimension of the global attractor [T], [5J, [TJ, 
|12j . [2H]. The Hausdorff and fractal dimensions of the global attractor satisfy the following 
estimate [5] (see also [6J, [39]): 

dim F ^ < Cl G 2/3 (log(l + G)) 1/3 , Cl = 01(7) 

which has been shown in [33] (following ideas of pQ) to be logarithmically sharp. 

If we accept the point of view that the small length scale can be defined as follows (see 

0, n, m, m) 



• (L2) 
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then up to logarithmic correction we have 

^^W- (1-3) 

This heuristic estimate for the small length scale is probably the best one can hope for 
since it matches, up to logarithmic term, the physically asserted estimates for the enstrophy 
dissipation length scale [30j . We also observe that the estimate (11.31) is extensive, that is, 
independent of the size of the spatial domain provided that its shape is fixed. 

Another rigorous definition of the small length scale can be given in terms of the number 
of determining modes, nodes, or volume elements (see |12| . |14j . |lbj . [25] and the references 
therein). It was shown that if N is sufficiently large and N equal squares of size ld n tile the 
periodic spatial domain, then any collection of points (one in each square) are determining 
for the long time dynamics of the 2D Navier-Stokes system. The best to date estimate 
for iV was obtained in [28J: 

N<c 2 G, 

where c 2 = 02(7) depends only on the aspect ratio 7 < 1. (An explicit estimate for c 2 was 
obtained in [26]: c 2 {i) = (68/(7vr)) 1 / 2 .) 

Therefore the small length scale defined in terms of the lattice of determining nodes 
satisfies 

^> C2 (7)- 1/2 ^- (1-4) 

We observe that this estimate is not extensive, that is, ld n scales like A _1//2 if VL is replaced 
by XQ, A > 0. 

We point out here that for the 2D Navier-Stokes system with analytic forcing the results 
of [IT], [18] provide the existence of a finite number iV of instantaneously determining 
nodes comparable with the fractal dimension of the attractor. These nodes, however, can 
be chosen arbitrarily (up to a subset of fl N with 2iV-dimensional Lebesgue measure zero) 
and therefore do not naturally define a regular lattice of determining nodes. 

The best to date estimate for the analyticity radius of the solutions of the Navier-Stokes 
equations with analytic forcing term / was obtained in 



/a - C3(7 W(l + ^G)V4- (L5) 

Relating the radius of analyticity to the dissipative small length scale (see also [23] in this 
regard) we note that up to a logarithmic correction the estimate (11.51) coincides with (11. 4ft . 
but both are worse than ( 11.211 . where the latter coincides, as we have already pointed out, 
with the physically asserted estimate of [30] . 

In this paper we focus on the 2D space periodic Navier-Stokes system with damping 

2 

d t u + u l diU = —fiu + vA u — V p + f, 

i=l 

divu = 0. 



;i.6) 



By adding the Coriolis forcing term to (II. 6p one obtains the well-known Stommel-Charney 
barotropic model of ocean circulation [I], [10], [35], [37] . Here the damping \iu represents 
the Rayleigh friction term and / is the wind stress. For an analytical study of this system 
see, for instance [5], [22], [23], [4"T] . and the references therein. In a follow up work we will 
be studying the effect of adding rotation (Coriolis parameter) on the size of small scales 
and the complexity of the dynamics of (11. 6p . Therefore, we will focus in this work on the 
system (II. 6p . We also point out that in this geophysical context the viscosity plays a much 
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smaller role in the mechanism of dissipating energy than the Rayleigh friction. That is 
why in this work the friction coefficient fi > will be fixed and we consider the system at 
the limit when v — » + . 

Sharp estimates (as v — > 0) for the Hausdorff and the fractal dimensions of the global 
attractor of the system (11 .6p were first obtained in the case of the square-shaped domain 
in [25] (7 = 1). Then the case of an elongated domain was studied in [27] (7 — > 0), where 
it was shown that 

dim F A <c 4 D, D — H rot /H°°M (1.7) 

where c 4 is an absolute constant (c 4 < 12). This estimate is sharp as both v — > and 
7^0. Therefore the small length scale defined as in (11.21) is of the order of 

ini \ 1/2 f iw \ 1/2 Ifil 1 / 2 



l *~\tim F A) ~ VII rot /IU J ~ DW ' (L8) 

This heuristic estimate is, in fact, a rigorous bound for the small length scale expressed in 
terms of the number of determining modes and nodes [26] : 

Idn = C5 (^] 1 = C 5 (jr-TTT-) 7 ' C 5 = 6§1/4 - ( L9 ) 



DJ ° VII rot /|| ( 

This means that any lattice of points in Q at a typical distance / < Idn is determining. 

The main result of this paper is in showing that the analyticity radius l a of the solutions 
of the damped-driven Navier-Stokes system (11.61) lying on the global attractor is bounded 
from below and satisfies the estimate: 

cm 1 / 2 

a ~ £>V2(1 + log£)l/2' (1 ' 10) 

which up to a logarithmic correction agrees both with the smallest scale estimate (11.81) 
and the rigorously defined typical distance between the determining nodes (II. 9p . 

It is worth mentioning that this point of view of relating the radius of analyticity of 
solutions on the Navier-Stokes equations to small scales in turbulence was also presented 
in [23]. 

This paper is organized as follows. In section [2] we employ the Gevrey-Hilbert space 
technique of [15J to derive a lower bound for the radius of analyticity of the order 

This bound considerably improves, for a fixed fi > 0, the lower bound (11.11) for the clas- 
sical Navier-Stokes system as v —>■ + (see also Remark 12.11) . Let us remark that as an 
alternative to the Gevrey regularity technique for estimating small scales one can apply 
the ladder estimates approach presented in [9] to obtain estimates for the small scales 
in (TIISD (see also [H]). 

In section E] the estimate (11.101) is proved for the system (11.61) following 



2. Gevrey regularity of the damped Navier-Stokes system 

As usual (see, for instance, [I], [6], [32], [38]), we write (II. 6p as an evolution equation in 
the Hilbert space H which is the closed subspace of solenoidal vectors in (L2(f2)) 2 with 
zero average over the torus Q = [0, L/j] x [0, L\. 

d t u + B(u, u) + vAu + fiu = /, u(0) = u . (2.1) 
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Here A = —PA is the Stokes operator with eigenvalues < Ai < A 2 < . . . , B(u,v) = 
P(Yli=i u l div) is the nonlinear term, / = Pf G H, and P : (L 2 (fi)) 2 — > H. 

We restrict ourselves to the case 7 = 1 and, in addition, assume that VL = [0, 2n] 2 (this 
simplifies the Fourier series below). The case of the square-shaped domain Q = [0, L] 2 
reduces to this case by scaling. Furthermore, any domain with aspect ratio 7 < 1 can be 
treated in the similar way, the absolute dimensionless constants ci, C2, ■ ■ ■ below will then 
depend on 7, however. 

A vector field u G H has the Fourier series expansion 



u 



Uj€^' X , Uj G C 2 , U-j = Uj, Uj ■ j = 0, u = 0, 

and 



3& 2 



M| 2 = NIL = E 



|2 



The eigenvalues of the Stokes operator A are the numbers |j| 2 , and the domain of its 
powers is the set of vector functions u such that 

(2tt) 2 ^ |j| 4 >/ = \\A a u\\ 2 < 00. 

For r, s > we define the Gevrey space D(e rAS ) of functions u satisfying 

(2vr) 2 e 2T|i|2 >i| 2 = \\e rA3 u\\ 2 < 00. (2.2) 

We suppose that the forcing term / belongs to the Gevrey space of analytic functions 
/ G D{e^ Al,2 A 1 ' 2 ), so that (2tt) 2 ]T \j\ 2 e 2 ^\u 3 \ 2 = \\e^ Al/2 A 1 ' 2 f\\ < 00 (2.3) 

for some o\ > 0. We set 

tp(t) = mm(v\\^ 2 t, ai). 

The norm and the scalar product in D(e <p ^ >Al/2 ) are denoted by || • || v and (-, respec- 
tively. 

We assume that uq G D(A 1 ^ 2 ) and take the scalar product of (12.11) and Au in D(e lp ^ Al/2 ) 

1 /2 

for sufficiently small t < 01/ (vA/ ). Since 

{e^ Al/2 d t u(t),e^ Al/2 u(t)) = l -d t \\A 1/2 <t)\\l - uX\ /2 (Au(t), A^u®)*, 
we obtain 

±d t \\A 1/2 u\\l + vWAuWl + ^WA^uWl = (B(u, u)Au) v + u\\ /2 (Au, A 1 ' 2 )^ + {A 1/2 f, A x l 2 u) ip . 

(2.4) 

Next we use the key estimate (see [15], [12], [ID]) for the nonlinear term in Gevrey spaces 



(\\Ai 
1 + l0g AX4i 



Au\\i 



/2 «ll 2 
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and use Young's inequality for this estimate and for the last two terms in (12.41) : 
d t \\A l '*u\\l+v\\Au\\l < 

WM\l \ , n , „, 1/2 „ a , U l/2 f\\% < 



<^P 1/2 ^li(l + logJg^) +2^11^,1^ 
<^||^||j(l + ]og- 



2/2 



\\M\l V , II^VIIJ 

-+- V Ai -+- 



.\\A l/2 »\\y 2a* ' 

where c 2 = 2c\ + 1. Next, using the inequality —az + (3(1 + log z) < (3\og(3/a (see 
we find 



-v\\Au\\l + | ll^l^l + log < f ||^||tlo g 



c 2 pV 2 M || 2 



and obtain the differential inequality 



■'-«||^ S — \\A-'-U\\ v lOg 

Hence the function 



dt\\A^uf v < ^il^H^log p» + + ^1_LN, 

v V v A\V l 2/i 



where lne = 1, satisfies 

d t y(t) < v\\Cj,y 2 logy, c 3 = max(l, c 2 /2). 
Therefore ?/(£) < 2y(0) for as long as 

t< (2u\ 1 c 3 y(0)\og2y(0))- 1 . 

In other words, 

WA^uWl < 2p 1 / 2 n || 2 + c,(u^f/ 2 \\A l l 2 fU + c,\ x v\ c 4 = e/c 2 , 
as long as < t < T*(\\A 1/2 u \\), where 

T 7 = <?r ,-\ ^ ^ll A1/2 ^ll 2 i ll A1/2 /ll^i | r \ w ^9 ^ll^ 1/2 "oll 2 i ll A1/2 /lk | ~ 
2c 3 z/Ai — ^ + Ai!/3 /2 M i/2 + ej log (^2 — ^ + Ail/3 / V / 2 + e 

We now observe (see Lemma 13 . 1 j) that on the global attractor or in the absorbing ball 
we have, respectively, 

\\A l ' 2 u(t)\\ < W Al/2f \ t e R, p 1/2 u(t)|| < 2 l|Al/2/l1 , t > ToiWA^uoW). 
\i // 

Therefore we have the following lower bound for T*\ 
In the limit — > + we have 

p i/ 2/ ||2 pl/ 2/ | U 



and we can write the lower bound for T* as follows 

T* > c 6 [z/Ai/^logD]" 1 
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where 



1^/11 II rot/Hi^ 1 wtae M/ | Un | 



In terms of the analyticity radius l a the lower bound for T* takes the form 

°- D 2 logB' 
Thus, we have proved the following theorem. 

Theorem 2.1. Suppose that f £ D(A 1 / 2 e CT1 " 4l/2 ) /or some o"i > 0. T/ien a solution u lying 
on the global attractor A is analytic with analyticity radius 

Iftl 1 / 2 



where 



Moreover, 



L ~ {(D* + D 1 + l)]og(D* + D 1 + l) ' ai 

D |rot/yg[ Di _ ||avvi| c 



c 8 



ftl 1 / 2 



/ a >^-L_ as _ +. (2.5) 
The constants c-j and c§ depend only on the aspect ratio of the periodic domain fl. 

Remark 2.1. We observe that the estimate (12.51) for the system (jl.fip is of the order 
z/ -2 log(l/z/) as far as the dependence on v — > + is concerned, while the estimate f ll.il) 
for the classical Navier-Stokes system is, in this respect much larger; namely, is of the 
order z/~ 4 log(l/V). 

However, the estimate (12.51) is not sharp and will be improved in the next section. As has 
been demonstrated in [31] the Gevrey-Hilbert space technique does not always provide 
sharp estimates for the radius of analyticity. The mechanism explaining this has been 
reported in [31] by means of an explicitly solvable model equation. 



3. Sharper bounds 

In this section we obtain sharper lower bounds for the analyticity radius l a . This is 
achieved by combining the z/-independent estimate for the vorticity contained in the fol- 
lowing lemma and the L p -technique developed in [21], [31] for the uniform analyticity 
radius of the solutions of the Navier-Stokes equations. We observe that similar technique 
has been earlier established in [2] for studying the analyticity of the Euler equations. 

Applying the operator rot to (jl.fip we obtain the well-known scalar vorticity equation 

d t uj + u ■ Vu = vAuj — jiu + F, (3.1) 

where ui = rot u, F — rot f,u = V^A" 1 ^, so that u ■ Vlj = V ± A~ 1 lj ■ Vlj = J(A~ 1 u, lj), 
and V x = (-St, , oi ), J(a, b) = V ± a ■ V6. 

Lemma 3.1. (See [26J.) The solutions u(t) lying on the global attractor A satisfy the 
following bound: 

M*)lk»< l|rot/lkfc , tern, (3.2) 

where 1 < k < oo. 
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Proof. We use the vorticity equation (13. II) and take the scalar product with ui 2k ~ l , where 
k > 1 is integer, and use the identity 

{J (if), <p), ^> 2k - 1 ) = (2k)- 1 J J (if), tp 2k )dx = (2k)- 1 J &v( V 2k V L ip)dx = 0. 

We obtain 

M\l*;%ML 2h + (2k-l)p J \\7tu\ 2 LU 2k - 2 dx + fi\\Lj\\i k 2k = 
= (rot/,^- 1 )< Hrot/iUMI^ 1 . 

Hence, by Gronwall's inequality 

\Ht)\\L 2k < Ik^ll^e-^ + ^Hrot/ll^a-e-^), 
and passing to the limit as k — > oo we find 

||w(t)||oo < ||^(0)|| oo e^ + /i- 1 ||rot/|| oo (l- e ^). 
Now, we let t — > oo in the above inequalities and obtain 

limsup||^(t)|| i2fc < l|rQt/l|Lafc , 1 < k < oo, 

which gives (13.21) since the solutions lying on the attractor are bounded for tel. □ 

As before we consider the square-shaped domain Q = [0,L] 2 and it is now convenient 
to write (11. 6p in dimensionless form. We introduce dimensionless variables x' , t' , u 1 and p' 
by setting 

x = Lx , t=(L 2 /v)t', u = (v/L)u, p=(u 2 /L 2 )p, /i = (z//L 2 )//. 
We obtain 



dsn' + u ' ld 'i u ' = ~^' u + A ' u ' - V V + /'> 



i=i 

divV = 0, 



(3.3) 



where x' G tt' — [0, l] 2 , /' = (L 3 /u 2 )f. Accordingly, the dimensionless form of (13. ip is as 
follows (we omit the primes): 

dtOJ + u ■ Wuj = Au — jjoo + F. (3.4) 

Remark 3.1. For dimensionless variables v! and &' the estimate (13. 2p with k = oo takes 
the form 

ll 'll ll +' Ml ^ ll rQt 71loo l|rot/||ooL 2 

oo = rot u oo < = = D. (3.5) 

/i' Z//i 

The next lemma is similar to the main estimate for the space analyticity radius in [31J. 

Lemma 3.2. Suppose that F is a restriction to Q (that is, y = 0) of a bounded x-periodic 
analytic function F(x + iy) + iG(x + iy) in the region \y\ < 5p and 

M 2 F = sup (F(x + iy) 2 + G(x + iy) 2 ). (3.6) 

xGQ, \y\<Sp 



Let p > 3/2 and let 



CM*/fi ■ 
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Here [and throughout) C is a sufficiently large universal constant and M 2p > ||^o||l2 P - 
Then the solution u(t) is analytic fort > and for < t < t the space analyticity radius 
of u(t) is greater than 

. .'t l l 2 1 1 1 

o[t) = mm -— , 



C ' Cpt( 2 P- 3 )/ 4 PM 2p ' Cpt( 2 P- 3 )/( 4 P+ 6 )M 2 2p/(2p+3) ' pt 1 ' 2 M 2p } 

Proof. We solve ( 13.41) by a sequence of approximating solutions (see [29], [31]). We set 
= and cu^ = 0. Then for u^ n \ u^ n > we have the equation 

<W n) - Aw (n) + w (n ~ 1} • W n) + fiuj (n) = F 

(3 7) 

£>>(0) = uj q = rot M , M (n) = V^-^W. 

The solutions tu^ and for t > have analytic extensions uj^ +i9^ n ' and u^ 1 +iv^ 
and since the system (13.7)) is linear, their analyticity radius is at least 5p. They satisfy 
the equation 

d t {uj {n) +i9 {n) ) - A(o; (n) +i9 {n) ) + (u (n_1) +iv {n ~ 1) ) ■ V(cu (n) +i0 {n) ) +fi{J n) +i6^ n) ) = F+iG, 
or, equivalently, the system 

<W n) - Acu (n) + uuj {n) + u (n_1) • W n) - • W (n) = F, 

(3 8) 

d t 9 {n) - M {n) + /ifl (n) + u^- 1] ■ W (n) + v^-V ■ Vcj (n) = G, 

where, as before, = V ± A^ 1 cu^ n \ v^ 1 ' = V^A -1 ^™- 1 , and the differential operators are 
taken with respect to x. In view of the analyticity of the solutions we have the Cauchy- 
Riemann equations 



dyj dxj 

&>) QQ{n) 

3 = 1,2, 



(3.9) 



dxj dyj 

and the similar equations for and v^ n \ 
Let e > 0. We consider the functional 

^ n (t)= f [ (u {n) (x,ats,t) 2 + 6 {n) (x,ats,t) 2 + e) p dxds. (3.10) 
Jo Jn 

We also set 

Q n (x, s, t) = uj {n) (x, acts, t) 2 + 6 {n) (x, acts, t) 2 + e. 

Here t G M + and a G M 2 . The combination ats will play the role of the variable y; p > 3/2, 
and £ > is arbitrary. 

We differentiate if) n (t) taking into account (13. 8p and use the Cauchy-Riemann equa- 
tions (13. 9p to handle the derivatives with respect to y. We obtain 

^-d t ^ n {t) + I = h + I 2 + I 3 + h, (3.11) 
2p 



where 



J = If Q^ 1 (|Va; (n) | 2 + |Ve (n) | 2 + //(a;(" ) ) 2 + A t(^" ) ) 2 )dxds+ 
Jo Jn 

+2(p-l) / / gr 2 (^ (n) V^ (n) +fl (n) W (n) ) 2 dxds, 
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and 

I 1= f / QP _1 (-w (n) V0 (n) + # (n W n) ) • asWs, 
Jo Jn 

[ [ Q p - l (uj^Vuj^ + 0^W (n) ) • u^dards, 
Jo Jn 

Jo Jn 

Jo Jn 

The arguments of Q„ are x, s,t, and the arguments of u^ n \ 9^ n \ uS n \ and are x, ats, 
and t. 

For an arbitrary rj > we have 



/ 2 
4 



h<r) f I Q p ~WV^ n) \ 2 + |W (n) | 2 )cfeds+ 
Jo Jn 

f [ QS" 1 (M n) ) 2 + (^ (n) ) 2 )l«rs 2 ^<^o + C',|a|V n (t). 
Jo Jn 



(3.12) 



Next, 



For i~3 we have 



J 2 = 2. / / VQ£ • u^'^dxds = 0. (3.13) 
2p Jo Jn 

h < Vh + C V f [ Ql \v^\ 2 dxds < r)I + (3.14) 
Jo Jn 

where 

/ pi p \ (p-i)/p 2 / f 1 f \ 1 ^ p 

I' 3 =l Q n (x,s,t) p2/ ( p - l) dxds) ? / 3=^(/ / |^ n_1) (a:,ats,t)| 2p da:ds) . 

\Jo Jn / „•_■, \Jo Jn / 



We write I' 3 as follows 



I's = \\Q P n /2 \\l g{ n o) , fi o = n x [0, 1] CM 3 , /? = 2p/(p - 1), 2 < < 6, 
and use in fi the Gagliardo-Nirenberg inequality 

\\A\\ L[jm < C7||^||K ) / 1V X> .^||^ + 6711^11^) 
for A = s) = Q?/ 2 (x, s, t). We have 
||V x ,s^4||i 2 (n ) = ||V x ,s Qn /2 |ll 2 (no) — 

P 2 [ [ Qr 2 ((^ (n) Vw (n) + d^Vd^) 2 + t 2 (^a • Vu (n) - u (n) a ■ V9^) 2 )dxds < 
Jo Jn 

< c P 2 (i + M 2 t 2 )/ . 



Hence, 

Next, \\Q p n /2 \\ 2 L2{ni)) = Mt) 



IIV.^H^ = \\V x , s Q^ 2 \\l% o) < C(l + \a\ 2 t 2 f^ir p . 



= ll« /2 |lSnor = ^) (2P " 3)/4P 
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and 
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I3 < C(i + \a\H 2 f 2p f /2p Mt) {2p ~ 3)/2p + CMt)- 



(3.15) 



We now consider 7g. Since vf" (x, 0,t) = (the solution restricted to y = is real- 
valued), we have (using the Cauchy-Riemann equations for Vj) 

2 ,1 

\vj (x, ats, t)\ 



f 1 

^^ctfcts / dy k Vj n (x, atsr, t)dr 
k=i Jo 

2 rl 

V^afcts / dkU ( j l ~ 1 \x,atsT,t)d,T 
k=i Jo 



Then 



2p \ 1 /P 

dxds I < 



7g = I / / V^afets / <9fcuj n 1 ^(x, aisr, £)dr 

i=i y° Jn fc =i ^° 

C\a\ 2 t 2 (J^ J J \Vu (n - 1 \x,atsT,t)\ 2p s 2p drdxd^j 

C\a\ 2 t 2 ( [ s 2p ds [ dr [ \Vu {n - l \x,atsT,t)\ 2p dx) 
\Jo Jo Jn J 



i/p 



Since u = V A uj, we have (see 120 



\Vu(x)\ 2p d3 



l/2p 



Therefore 



I3 <Cp 2 \a\H 2 



-1 r l 




-'p 



\J n ' l) (x, octsr, t)\ P dx s 2p dsdr ) < 



Jo Jn 



Cp 2 \a\ 2 t 2 (£ J Ql^dxd^j ' <Cp 2 \a\ 2 t 2 ip n - 1 {t) l/p , 

where we have used h(sT)s 2p dsdT < (2p)^ 1 h(s)ds. Combining this with (13.141) 

and (13.151) we obtain 



h < 



V % + C»|V(1 + |a| 2 t 2 ) 3 ^7 3/2 ^n(t) (2p - 3)/2 ^n-l(t) 1/P + CpMH^itfl^it) < 

V I + C v p 2 (\a\t)^- 3 \1 + |a| 2 t 2 ) 3 /^- 3 )^n-i(t) 2/(2 ^ 3 Vn(t) + Cp 2 \a\H 2 ^M /P Mt). 

(3.16) 



Finally, we estimate I4. 



I A < 




JO. 



QP- 1 ^) 2 + 6^) 2 ) w + (F 2 + G 2 )/(4 w ))dxds < 



V I + C v (M 2 F /v) 
where Mp is defined in (13. 61) . 




(3.17) 



Ql~ l dxds < V I + C v (M 2 /fi)Mt) {p ' 1)/p , 



JO, 
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Taking r) > sufficiently small we infer from f[3TTTD . f[3TT2]) . fl37L3l) . f[3T6D . f[3~T7j) 

d t Mt) < Cp\a\ 2 Mt) + Cp 3 |a| 4 ^^- 3 )^/^- 3 )^-i(t) 2/(2p - 3 Vn(t) + 

C7p 3 1 a | (#+6)/(2 P -3) t (4 P+ 6)/(2 P -3) ^ (f )2 /(2p-3) ^ (f ) + 

Cp 3 |a|VVn-i(0 1/p ^) + Cp^t^-^Ml/ii, 

where 

^„(0) = / (cu (x) 2 + e)^x. 
./n 

We set (f n (t) = ip n (t) 1 ^ 2p and obtain the differential inequality for cp n : 

d t <p n (t) < C7|a|V»(*) +Cp 2 | a |^/^- 3 )^/^- 3 )^_ 1 (0 4p/(2p - 3 Vn(*)+ 

Cp 2| a |(4 P+ 6)/(2 P -3) t (4 P+ 6)/(2 P -3)^_ i(t) 4 P /(2 P -3)^ ri(t) + 

We now use the Gronwall-type Lemma [3.31 from [3T] below and see that tp(t) < 2ip(0) 
on the time interval specified in (13.181) . (13.191) . and letting £->0we obtain 

-1 p 

(J n \x,ats,t) 2 + 9 {n) {x,ats,t) 2 ) P dxds < 2 2p Mj, 



'0 Jn 

for t > 0, lali < <5f and 



t < min(Ai, A 2 , A 3 , A 4 , A 



5), 



where 



Ax 

^3 

A, 
A 5 



C\a\ 



1 



Cp(4p-6)/(6 P -3) | tt |4 P /(6 P -3) M2 4 PA 6 P- 3 ) ' 
1 

C p (4 P -6)/(6 P +3) | a | (4 P +6)/(6 P +3)jVf 2 4 P/( 6p+3 ) 
1 



^2/3| a |2/3 M2 2 / 3 



CM 2 F fi~ l 



We now set 



Then the condition 



CM 2 F /fi 



t < min(A u A 2 , A 3 , A A ) 
can be written in terms of y = at as follows 

.'/' 2 1 1 

\y\ < min 



C ' Cpt^-^l A PM 2 p Cpt( 2 P- 3 )/( 4 P+ 6 )Mj /(2p+3) ' pt 1/2 M 2p 
Now for to defined in (I3.20p and 

't 1 ' 2 1 1 1 



5{t) = min 



C ' Cpt^V^M^ Cpt( 2 P-m*p+G)M% /i2p+3) ' pt 1 / 2 M 2p ' 



6 r 



(3.18) 



(3.19) 



(3.20) 



(3.21) 



(3.22) 
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we have for < t < t and \y\ < 5(t) 



(J n \x,sy,t) 2 + 9 (n) (x,sy,t) 2 )dxds < 2 2p M 2p 



o Jn 

2 



for all integer n > 1. Therefore for any y e R with |y| = 1 this gives that 
,<5(t) /• 

/ (cu (n) (x, sy,t) 2 + # (n) (x, sy,t) 2 )dxds < 2 2p 5(t)M 2p 
o in 

and since f(sy)ds < B, \y\ = 1 implies J^ <s f(y)dy < 2rr5B, we obtain 
/ / (uj^\x,y,t) 2 + 9 {n) {x,y,t) 2 )dxds <2n2 2p '8 (t) 2 M 2p . 

J\y\<8(t) Jn 

This estimate is uniform in n and as in |21j . [31] we obtain the existence of an analytic 
solution of (13.41) with analyticity radius satisfying (I3.22[) . The proof is complete. □ 

Lemma 3.3. (See [31J.) Let y n (t) 6 C 1 [0,T] be a sequence of non-negative functions 
satisfying yo(t) < M for < t < T, and y n (0) < M for n > 1. Suppose that on the 
interval < t < T 

N 

d t y n (t) < J2 K J ta3 yn(tpy n -i(ty>, 

i=i 

where Kj > 0, ctj > —1, (3j G R, and 7j > are given constants. Then y n {t) < 2M for 
all n = 0, 1, 2, . . . provided that 

/ / a,- + 1 
< t < mm I T, mm t 

j=i,...,N y NKfiPi +^Mft+^- 1 

where (3 + = max(/3, 0). 

We can now state the main result of this section. 

Theorem 3.1. The solutions on the 2D space-periodic damped- driven Navier-Stokes sys- 
tem (1 1 . /?/mg on t/ie global attactor A are analytic with space analyticity radius l a satis- 
fying the lower bound 

CD*{1 + V*D)W D -^f^- (3 - 23) 

Proof. We first observe that (13.231) is equivalent to the estimate 

k ~ 671)1/2(1+ log (3 - 24) 

for the equation written in dimensionless form. 
Next, by Young's inequality 

pt (2 P -3)/4p M < CpM 4p/(A P -3) t l/2 + 
pt (2 P -3)/(4 P +6) M 2 P /(2 P +3) < C pMt 1/2 + t" 1/2 . 

Hence, the estimate (13.221) can be written as follows 

/ 1 1 / 2 1 \ 
<$(*)> min — , . ... ,5 F . (3.25) 
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The solutions lying on the attractor are bounded in L 2p : 

K*)||j^<Map, M 2p <CM oc . 

Setting 

p = C(l + \ogM 00 ) 

we see that 

p ( M J/(4 P -3) + < + logMoo)Moo 



<2p 

and therefore 



p=C(l+logMoo) 



' t l/2 

6(t) > min ( — - 



C ' C(l + log 



At the moment of time 

1 



C(l + logMoo)M 



which for sufficiently large (the case of our interest) is smaller than t defined in (13.201) 
(the details are given below) we have 

5(f) > 1 



CM^l + logMoo) 1 ^ 
Since < D (see ( 13.51) ). it follows that 

1 



6(f) > 



CD 1 /2(l + logD)l/2 



By the invariance property of the attractor we see that on the attractor the above estimate 
holds for all £*, which proves (I3.24p . 

To complete the proof it remains to show that 

** < ^ = , (3.26) 



C(l + log D)D ~ u CM|,/a« 

where in the expression for to we reverted to the prime notation for the dimensionless 
damping coefficient // and the forcing F' . We relate the forcing term and its analytic 
extension by the equality 

M F , = iqrot'/IU, K = K(F,6 F ). 

Recalling that /' = (L 3 /z/ 2 )/, /jf = (L 2 /u)fi, and x' = (1/L)x we see that 

v 

to = CK 2 /iL 2 ' 
Hence, f !3.26|) goes over to the condition 



C(l + \ogD) > 



which is obviously satisfied for all sufficiently small v > 0. The proof is complete. □ 
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4. Concluding remarks 

We have shown that the solutions lying on the attractor of the 2D space-periodic 
damped-driven Navier-Stokes system, the Stommel-Charney barotropic model of ocean 
circulation without rotation, with analytic forcing have space analyticity radius which up 
to a logarithmic term coincides with the small scale estimates both in terms of the sharp 
bounds for the fractal dimension of the global attractor, and in terms of the spatial lat- 
tice of determining nodes. The derivation of this lower bound for the analyticity radius 
essentially uses the techniques developed in [3T] . 
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